Coherent control with shaped femtosecond laser pulses applied to ultracold molecules 
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We report on coherent control of excitation processes of translationally ultracold rubidium dimers 
in a magneto-optical trap by using shaped femtosecond laser pulses. Evolution strategies are applied 
in a feedback loop in order to optimize the photoexcitation of the Rb2 molecules, which subsequently 
undergo ionization or fragmentation. A superior performance of the resulting pulses compared to 
unshaped pulses of the same pulse energy is obtained by distributing the energy among specific 
spectral components. The demonstration of coherent control to ultracold ensembles opens a path 
to actively influence fundamental photo-induced processes in molecular quantum gases. 
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I. INTRODUCTION 

The interest to study interactions in ultracold molec- 
ular gases is propelling various research programs to 
produce and trap dense ultracold molecular ensembles 
IE El [J IE IE 0- The goals and perspectives in this 
field range from quantum computation y| and quantum 
scattering to ultrahigh precision spectroscopy 01 an d 
coherent ultracold chemistry The manipulation of 

molecular wavepacket dynamics with shaped femtosec- 
ond pulses may open exciting possibilities to this field. In 
recent years, the control of molecular processes by shaped 
fs-laser pulses through application of evolution strate- 
gies in a feedback loop has attained considerable success. 
Since it was proposed coherent control with shaped 
femtosecond laser pulses has been applied to a large va- 
riety of experiments, i.e. selective frag mentation |l3j . 
population transfer 14], ionization and high har- 
monic generation |16| . By using an iterative process, the 
method enables the discovery of an optimal laser pulse 
shape which drives the system towards a desired target 
state. It even became possible to decipher the molec- 
ular dynamics underlying the laser-molecule interaction 
by analyzing the optimal pulse shapes 0, . 

Theoretical work on photoassociation with pulsed 
lasers has started several years ago [l8( and first experi- 
ments with thermal mercury atoms in a gas cel l Il9ll and 
with sodium atoms in a magneto-optical trap [20( have 
been performed. More recently, optimal control calcu- 
lations were done proposing efficient ultracold molecule 
formation pH and vibrational cooling by appropri- 
ately shaped laser pulses. However, up to now there are 
no experimental results on the application of feedback- 
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controlled optimization to photoassociation or manipu- 
lation of ultracold molecules. In this article we present 
first experimental steps towards this goal by investigating 
the interaction of femtosecond laser pulses with ultracold 
ground state molecules in view of a pulsed photoassocia- 
tion experiment. The low temperature of the atomic and 
molecular sample, while not directly influencing the in- 
teraction with the femtosecond pulses, leads to molecules 
in very high lying vibrational levels and allows us to 
study molecules under the same conditions of typical pho- 
toassociation experiments. In this work we demonstrate 
in particular the enhancement of the electronic excita- 
tion and subsequent fragmentation of ultracold molecules 
from a magneto-optical trap using shaped femtosecond 
laser pulses, which were adaptively optimized by an evo- 
lutionary algorithm. 



II. EXPERIMENTAL SETUP 

In our experiment, schematically depicted in Fig. ^ 
about 10 7 85 Rb atoms are captured in a magneto-optical 
trap (MOT) at a density of 10 10 atoms/cm 3 and a tem- 
perature of 100 fi K. Diatomic rubidium molecules con- 
tinuously form in the magneto-optical trap due to either 
three body collisions or light assisted two-body collisions 
of trapped Rb atoms. They populate predominantly the 
highest vibrational states below the dissociation limit in 
the triplet electronic ground state [IE 0] • These Rb2 
molecules, which are not trapped by the trap magnetic 
field, are detected in this experiment via resonant two 
photon ionization and time-of-flight mass analysis. The 
ionization laser is a 15 Hz Nd:YAG pumped dye laser at 
wavelengths between 600 and 610 nm, a spectral width 
of 0.5 cm -1 and a pulse energy of 20 mJ. In the steady 
state of molecule formation and loss a maximum count 
rate of 0.5 Rbj" molecular ions per laser pulse is observed, 
measured at a wavelength of 602.5 nm. 
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FIG. 1: Experimental setup for iterative closed-loop maxi- 
mization of ultracold Rb2 excitation from the ground elec- 
tronic singlet or triplet states by shaped femtosecond laser 
pulses. The ultracold molecules are formed in a magneto- 
optically trapped gas of rubidium atoms. The laser pulse 
shapes represent individuals of an evolutionary algorithm. 
Their fitness is evaluated by the reduction of the RbJ sig- 
nal resulting from resonantly enhanced photoionization of the 
Rb2 molecules. 
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FIG. 2: Reduction of the RbJ molecular ion signal as a func- 
tion of energy of transform limited femtosecond pulses. The 
pulses have a central wavelength of 800 nm and 10 nm FWHM. 
The Dl atomic resonance at 795 nm is filtered out of the pulse. 
The ionization laser for REMPI is set to 602.5 nm. The solid 
curve represents a fit to a single photon excitation model as 
described in the text. 



The femtosecond (fs) laser pulses are generated in a 
TkSapphire oscillator (Tsunami; Spectra Physics) that 
provides pulses of 120 fs duration (FWHM) at a rate of 
80 MHz, a spectral width of AA = 10 nm and an energy 
of up to 18 nJ per pulse. Low pulse energies ensure that 
the laser-molecule interaction is well described in a per- 
turbative picture and no high intensity effects have to 
be considered. To modify the spectral components of 
the pulses we use a pulse shaper that allows indepen- 
dent phase and amplitude modulation [2^. It consists 
of a liquid crystal modulator (CRI; SLM-256) [H with 
2x128 pixels, placed in the Fourier plane of a double grat- 
ing zero dispersion compressor. A lens focuses the beam 
to a spot of 150 /jm diameter at the center of the trap, 
illuminating about 10% of the cloud volume. For the ex- 
periments the fs-laser is tuned in the range between 780 
and 820 nm. 



III. LASER-INDUCED LOSS 

Transform-limited fs-pulses with wavelengths near the 
Rb resonance lines and focused into the MOT are found 
to strongly interact with the trapped atoms, as observed 
through a significant decrease of MOT fluorescence. At 
small pulse energies, the trap loading rate outside the 
femtosecond beam can partly compensate for the losses 
in the focus, resulting in an area of reduced fluorescence 
predominantly where the beam is passing, whereas for 
high energies the trap is completely depleted. We at- 
tribute this to photon scattering from the femtosecond 
pulses, which causes an outward directed light force on 
the atoms, leading to trap loss. The effect is significantly 
stronger when the D2 atomic resonance at 780 nm was 



part of the pulse spectrum than for the Dl resonance at 
795 nm, reflecting the higher transition dipole moment 
of the D2 line. A narrow slit, representing a bandpass 
of 0.3 nm width, is scanned through the Fourier plane 
of the pulse shaper, which showed that only components 
resonant with the Dl and D2 atomic transitions are re- 
sponsible for this effect. Singly-charged rubidium ions 
have been detected when the central pulse frequency was 
resonant with one of the atomic transitions, indicating 
that resonant three-photon ionization contributes to the 
loss of atoms from the MOT. To study the laser pulse in- 
teraction with rubidium molecules, the atomic resonance 
components were removed from the pulse spectrum by a 
notch filter, realized by a physical block in the shaper's 
Fourier plane. In this way atomic losses from the MOT 
could be reduced below the detection threshold. 

The rubidium dimers interact with the fs-laser pulses 
over the entire accessible range of central wavelengths 
from 780 nm to 820 nm. As shown in Fig. [21 the molec- 
ular signal decreases rapidly at small pulse energies and 
levels off to 25% at a pulse energy of 0.6 nJ. As only 
molecules in the electronic ground state are detected, 
the signal reduction can be attributed to excitation by 
the fs-pulses. The process can be modeled by a simple 
rate equation for the number of detected ground state 
molecules: 

■ = Rfoim — (Pfs + rioss)A r Mol (1) 

where i?f or m is the production rate of molecules from 
trapped atoms, Tf s is the rate of molecular loss in- 
duced by interaction with the femtosecond laser, and 
Tioss the molecular loss rate induced by other processes 
such as rest gas collisions or gravitation accelerating the 
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molecules out of the detection volume. In the steady 
state dNyi \/dt — the number of molecules is given by 
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The curve in fig. represents a fit to the data assuming 
a linear dependence of T& on the pulse energy of the 
laser, i.e. T& oc Pf s . The assumption of a quadratic or 
higher order dependence of the excitation rate on pulse 
energy does not fit the data. This indicates that, in the 
regime of pulse energies employed in the experiment, the 
interaction with the molecules has the character of an 
effective one- photon excitation p7j . 

According to [23|. the molecules in the MOT initially 
populate the highest levels in the a 3 S+ state. Due to se- 
lection rules and Franck-Condon factors, they are prefer- 
ably excited to the 0~ and l g 5s5pi/ 2 states (see Fig. 
[3}. Emission back into the electronic ground state could 
only lead to signal reduction if the vibrational level pop- 
ulation moves out of the excitation window of the ioniza- 
tion laser. However, scans of the detection laser with and 
without fs-beam show reduced but qualitatively similar 
spectra which should not be the case for a vibrational 
redistribution. Instead it can be expected that excited 
molecules absorb further photons, so the whole process 
of molecular loss can be regarded as a resonance enhanced 
multi-photon excitation, followed by dissociation, predis- 
sociation or ionization. This happens either within one 
pulse, or, as the laser repetition rate is comparable to 
the lifetime of the first excited state, it occurs in the 
subsequent pulse. At high energies all molecules in the 
laser focus are excited or dissociated and the residual 
signal in Fig. [21 is due to molecules which did not inter- 
act with the femtosecond laser. This shows that most of 
the molecules are produced within a small volume inside 
the MOT which is consistent with the picture that they 
form at the MOT center where the atom number density 
is at its maximum [28| . At high energies all molecules in 
the laser focus are excited or dissociated and the residual 
signal in Fig. [21 is due to molecules which did not inter- 
act with the femtosecond laser, indicated by the dashed 
line. This shows that most of the molecules are produced 
within a small volume inside the MOT which is consis- 
tent with the picture that they form at the MOT center 
where the atom number density is at its maximum |28| . 
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FIG. 3: Potential curves of the rubidium dimer including spin- 
orbit interaction. Initially, the molecules are expected to pop- 
ulate the highest levels in the a 3 Ej [23l l24i . The molecular 
excitation by the femtosecond pulses is indicated by the ar- 
row. The shaded areas show the initial and final distribution 
of molecular vibrational states. 



algorithm. To reduce the search space for the learning 
algorithm we chose a mixed scheme of parametric am- 
plitude and free phase optimization where the algorithm 
tries to find the optimal pulse shape under the restriction 
that only a few sharp spectral peaks contribute to the 
pulse shape. During an optimization the parameters of 
these peaks, their spectral positions and amplitudes, are 
altered to find the most efficient excitation pulse. More- 
over, the phase was optimized freely in order to allow a 
temporal modulation of the pulse. The evolutionary al- 
gorithm administrates 31 individuals, each representing 
the pulse shaper parameters to produce pulses consisting 
of up to eleven Gaussian peaks of 7 cm -1 FWHM. 

The adaptive algorithm was applied to manipulate the 
excitation pulses with the aim to minimize the molecu- 
lar signal from the MOT. For each iteration the ion sig- 
nals corresponding to the best and worst individuals are 
protocoled together with the mean fitness of the whole 
generation. As depicted in Fig.^a), all three signals de- 
crease during the particulate optimization to about 70% 
of the initial value after 20 iterations. The spectra of 
the final best individuals of two successive runs shown in 
Fig. 0fb) display several peaks which coincide in some 
but not all spectral positions. The frequency span of the 
fs-pulsc supports our assignment of excitation to the 0~ 
IV. APPLICATION OF COHERENT CONTROL and l g 5s5pi /2 states (see Fig. EJ. By comparing the 

excitation yield of the best individuals with transform- 
limited pulses of the same energy it is observed that the 
optimized pulse excites the molecules on average 25% 
more efficiently, which demonstrates the feasibility and 
potential of adaptive control. 

We attribute the observed excitation enhancement to 
an increased spectral intensity at particular molecular 
resonances found by the evolutionary algorithm. Starting 
from a narrow band in the a 3 £+ state @[^|, molecules 
are excited into bound states below the Dl resonance. 



In order to demonstrate the practical applicability of 
coherent control concepts to ultracold molecules, the 
RbJ signal acts as an input for a self-learning optimiza- 
tion algorithm which autonomously programs the pulse 
shaper in a closed loop experiment. The algorithm is 
based on evolution strategies and is described in detail 
in (2^|. Because of the small molecular ion count rate 
and hence the low signal-to-noise ratio the signal is aver- 
aged over 128 dye laser pulses for each individual of the 
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FIG. 4: (a) Molecular ion signal resulting from the best, the 
worst and the mean individual of the population for each gen- 
eration during a closed loop experiment, (b) Femtosecond 
laser pulse spectrum of the final best individuals of two suc- 
cessive optimization runs under equal conditions with similar 
final optimization result. 



By shifting the peak positions, the algorithm finds tran- 
sition frequencies from this band to certain vibrational 
states, thereby sharing the pulse energy more efficiently 
than a broad Gaussian pulse. The algorithm therefore 
has a large number of possible solutions to choose from 
and so the final pulse shapes after an optimization are 
not identical, fn the spectral region between 12000 and 
12500 cm -1 , the vibrational level separation is about 10 
cm -1 in the 0~ and l g 5s5pi/2 states, respectively. The 
high density of states also explains the limited poten- 
tial of the optimization because the optimization factor 
depends on the chosen peakwidth which is limited by 
the shaper resolution. The Franck-Condon factors may 
also be relevant for the excitation process since they dif- 



fer for different vibronic transitions and favor particular 
frequencies which are enhanced in the experimentally ac- 
quired spectra. Yet, as the initial ground state popula- 
tion distribution in the vibrational states is not known 
accurately, no quantitative treatment or assignment can 
be made. 

V. CONCLUSIONS 

In this article we demonstrate the application of iter- 
ative adaptive control techniques to the manipulation of 
ultracold molecules. The minimization of molecular sig- 
nal from a rubidium magneto-optical trap results in pulse 
shapes that significantly excite more Rb2 molecules than 
transform-limited pulses of the same energy. The result- 
ing pulse spectra are not unambiguous because of the 
large number of possible optimal solutions within the ex- 
perimental accuracy. For future applications of shaped 
fs-pulses in photoassociation experiments, the dissocia- 
tive effect of the pulses has to be suppressed. On the one 
hand, photoassociated cold molecules may be dissociated 
by subsequent femtosecond pulses due to the high repe- 
tition rate of the femtosecond laser system. Therefore, 
photoassociated molecules have to be detected directly 
after their formation by increasing the repetition rate of 
the detection system, ideally matching the rate of the 
photoassociating beam. On the other hand, the excita- 
tion of ground state molecules may be coherently sup- 
pressed while enhancing the molecule formation through 
actively shaping the pulses. For this purpose, we cur- 
rently theoretically investigate realistic scenarios for fem- 
tosecond photoassociation of cold molecules by numeri- 
cally simulating the dynamical behavior of an atom pair 
interacting with an arbitrarily shaped electromagnetic 
field. It is also envisioned to use specifically shaped fem- 
tosecond laser pulses for the cooling of vibrational molec- 
ular excitations by pump-dump processes via an inter- 
mediate excited state. This approach could be applied 
to produce quantum-degenerate molecular gases in arbi- 
trary vibrational states or superpositions of states. 
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